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Abstract 

We investigate Kolmogorov wave turbulence in QCD or, in other words, we 
calculate the spectrum of gluons as a function of time, fk(t), in the presence 
of a source which feeds in energy density in the infrared region at a constant 
rate. We find an early, an intermediate and a late time form for the gluon 
spectrum. Wave turbulence in QCD turns out to be somewhat different than 
the turbulence in the case of -type theories studied by Zakharov, L'vov and 
Falkovich. The hope is that a good understanding of QCD wave turbulence 
might lead to a better understanding of the instability problem in the early 
stages of the evolution after a heavy ion collision. 
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1 Introduction 



One of the key questions in heavy ion physics is how rapidly and by what mechanism 
equilibrium is reached after a collision. Phenomenological analyses of experimental 
data at RHIC suggest a rapid thermalization [1, 2] but, so far, no convincing theoret- 
ical picture has emerged which naturally gives such a quick approach to equilibrium. 

The "bottom-up" picture of equilibration [3] is a detailed picture of how the 
initially produced hard gluons (at or near the saturation scale of the colliding ions) 
lose energy by radiating softer gluons. After they become sufficiently numerous 
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these softer gluons equilibrate amongst themselves and continue to absorb energy 
from the initially produced hard gluons. After the hard gluons have lost all their 
energy the system has reached complete equilibration. Although difficult to evalu- 
ate precisely, the time for this equilibration was estimated to be on the order of 3 
fm [5]. While it is an attractive picture the bottom-up scenario has a serious flaw 
in that the initially produced hard particles quickly acquire an instability as they 
spread out along the axis of the heavy ion collision as it was recently worked out by 
Arnold, Lenaghan and Moore [6] and was long ago advocated by Mrowczynski [7] 
(For other early disscusions, see Refs. [8-10]). This instability then becomes the 
dominant mechanism for the creation of softer gluons [6], more important than the 
Bethe-Heitler radiation used in the bottom-up picture [3,4]. The instability poten- 
tially speeds up equilibration, but so far it has been difficult to follow analytically 
how the QCD gluonic system evolves toward equilibration in the presence of the 
instability [11, 12]. Numerical studies [13-18] seem to indicate that the instability is 
effective at early times, but becomes less prominent in non-Abelian theories when 
the occupation number of low momentum modes becomes large. So far there is no 
analytic understanding of the numerical results. 

The phase space spectrum found in numerical studies by Arnold and Moore [13] 
has a resemblance to that of the Kolmogorov spectrum in turbulence [19]. Indeed, 
the problem of wave turbulence with an infrared source of energy discussed by 
Zakharov, L'vov and Falkovich (ZLF) [20] (for a nice overview, see [21]) would seem 
to have much in common with the instability problem currently under discussion for 
a non-expanding QCD medium and where the source for the instability is a collection 
of hard particles having a fixed asymmetry in their momentum distribution. Wave 
turbulence is a somewhat easier problem to deal with since one can take a spherically 
symmetric source of inflow of energy in low momentum modes. One can hope that 
a good understanding of QCD wave turbulence will lead to a better understanding 
of the early stages of evolution after a high-energy heavy ion collision. 

Wave turbulence in QCD appears to be somewhat different than the situation 
studied by Zakharov, L'vov and Falkovich [20]. The essence of the ZLF discussion 
is that waves, or particles, interact with each other locally in momentum. This is 
certainly the case, say, in a 4 theory where the scattering cross section for a high 
momentum particle on a low momentum particle is very small. Soft quanta in <p A 
theory interact very weakly with hard quanta because of the small overlap in their 
fields due to Lorentz contraction of the harder quanta. Soft gluons, on the other 
hand, interact with the currents of the harder gluons, and these currents do not 
decrease with the momentum of the harder gluons. In addition a gauge theory has 
a lower cutoff in frequency, the plasma frequencly, which has no analog in a theory 
of the 4 variety. Because of the strong interaction between soft and hard modes 
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along with the plasma frequency infrared cutoff QCD turbulence seems to be quite 
different than the ZLF problem. 

In the problems considered by ZLF a source, say, of energy inserts energy at an 
infrared scale and a sink extracts that energy at a higher scale. The flow of energy 
then proceeds from a low momentum scale, through intermediate energy scales and 
exits at some high energy scale. A closely related process omits the high energy sink 
and has the energy flowing from low to high momentum through all intermediate 
scales. What is different in the QCD case is that energy can be absorbed by very high 
energy particles from a low energy source without the necessity of passing through 
intermediate scales. Particle number, however, cannot be directly transferred from 
low to high energy but must flow through intermediate scales, or be created at 
high energy through inelastic reactions. Thus in QCD there are important, even 
dominant, interactions which are not local in momentum. 

In this paper we calculate the analytic form of the spectrum of gluons in the 
presence of a source which feeds in energy density at a constant rate e = m\ja. 
We suppose the energy comes into our system uniformly in space in the form of 
gluons which have a spherically symmetric momentum distribution and which are 
inserted uniformly in phase space just above the plasma frequency cutoff in a range 
m < uj < m with m the plasma frequency and m on the order of m. mo is the single 
dimensionful parameter in the discussion given in Sec. 2 while in Sec. 3 we allow 
the incoming energy to be spread uniformly in phase space in a region < k < k 
of momenta in which case k is a separate dimensional parameter if we choose k 
to be a scale larger than and independent of m . We always suppose a, the gluonic 
coupling, to be small. 

At late times after the source has become active, m t > (I/a) 9 / 5 in case the 
source energy is deposited in m < uj < m, the system of gluons is very close to 
thermal equilibrium. The incoming energy is transferred from the scale m to the 
scale given by the temperature T by direct absorption of soft gluons in an inelastic 
3^2 process and, parametrically equally as important, by elastic scattering of soft 
gluons (scale m) on hard gluons (scale T) with the soft gluons losing energy to the 
hard gluons. Both m and T are slowly increasing with time. 

In the time domain (I/a) 7 / 5 < m t < (1/a) 9 / 5 gluons having momentum much 
greater than m are in thermal equilibrium while gluon occupation numbers are 
significantly above the thermal curve in the domain m < uj < m. 

Finally, in the early time domain, 1 <C m^t < (1/a) 7 / 5 , the system is far from 
thermal equilibrium in both the high and low momentum regimes. If po(t) is the 
maximum scale to which gluons have evolved, then f po ^> 1 while = in the 
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domain m < w < This latter distribution is a sort of equilibrium distribution 
although it does not match onto a genuine equilibrium distribution at the scale p . 

While intermediate momentum scales obey ~ 1/k in all time domains, so long 
as m t ^> 1, there is, nevertheless, a flow of energy and particle number from soft 
to hard scales through these intermediate scales. In the case of energy this flow is 
negligible compared to the direct transfer of energy from soft to hard scales, while 
in the case of particle number flow it is one of the dominant mechanisms of gluon 
number growth at the hard scale. 

The dynamics which we use throughout this paper is the Boltzmann equation 
with a collision term consisting of 2 <-> 2 and 2^3 gluon processes. We have 
assumed the absence of important longrange coherent fields which, if present, could 
create problems for our approach. 



2 The spectrum of gluons as a function of time 



The problem we are concerned with here is very similar to the problem of wave 
turbulence discussed by Zakharov, L'vov and Falkovich [20]. However, Zakharov et 
al. consider scalar theories where interactions among particles are relatively local 
in momentum while the case of interest here, QCD, allows important interaction 
between low momentum and high momentum particles. In addition in a gluonic 
theory there is a minimum allowed frequency, the plasma frequency, at which waves 
can propagate while a corresponding phenomenon does not exist with purely scalar 
interactions. 

At time, t, equal to zero we turn on a source of gluons which feeds a rate of 
increase of energy density e = = which is constant in time and uniform 
in space. The energy is incoming isotropically in low momentum modes but as the 
density of gluons increases in time the source is modified so that the incoming gluons 
always have energy just above the plasma frequency, say m < u < m where m is 
the plasma frequency and m is of the same order as m. The question we wish to 
answer concerns the occupaton number, of gluons as a function of time. As we 
shall see, the form of the spectrum changes over the course of time having an early 
time, an intermediate time and a late time form. The late time form of is the 
simplest and the most straightforward to describe so we begin there. Throughout 
we limit our discussion to the case of a fixed QCD coupling, a. 
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2.1 The spectrum for m t > (l/a) 9/5 



At sufhcently large time it is apparent that the system will be extremely close 
to thermal equilibrium since the source can only be a very small perturbation on 
a system having a lot of energy. (The scaling solution for constant energy flow 
found by ZLF [20] (see also [21]), fk ~ (m//c)~ 5 / 3 , connects an infrared source of 
energy to a quantum equilibrium solution of the Boltzmann equation which exists 
for k > ko whwere fk Q — 1. The non equilibrium solution in the intermediate region 
of k is necessary in a 4 -type theory considered by ZLF where energy flow is local in 
momentum. In the QCD case energy is transferred directly from the low momentum 
source to the high momentum modes allowing the intermediate momenta also to 
be in equilibrium. Indeed, we shall see throughout this paper that ~ m/k is 
required, in intemediate momentum regions, in order that the rate of emission and 
absorbtion of gluons of momentum k by harder gluons exactly cancel.) Thus energy 
conservation gives 



rrin 



e(t) = ^ • t = g E T 4 (1) 
a 



with g E = 2(JV= - 1)%, or 



r = m [^V /4 . (2) 



g E aJ 



The corresponding plasma frequency is 



An 

uj P = m = y —aN c T. (3) 

and is related to the Debye mass by m = ttid/V^- At large times the energy 
source inserts gluons into the system just above an energy m which is far below the 
temperature of the system T. What we now wish to describe is how the energy gets 
transferred from the scale m to the scale T which is where almost all the energy is 
located. First we remark that, for very small a, the spectrum is given by 

which, when cu/T <^ 1, can be approximated as 

h * -■ (5) 

When fc/m > 1 we shall not distinguish between u and k in which case (5) could 
as well be written as 

h - \. (6) 
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Of course if the system were in exact equilibrium there would be no transfer of 
energy and (6) would represent a system in complete equilibrium without flow of en- 
ergy or particle number. Because we have a source of incoming energy the occupation 
number will be slightly higher than the equilibrium solution in the low momentum 
region. As we shall now demonstrate these low momentum gluons directly transfer 
energy to the hard gluons, having momentum T, with only a small fraction of energy 
flowing from the scale m through an intermediate scale k,m <C k <C T, and then to 
the scale T. However, particle number, which must increase as T 3 , cannot be trans- 
ferred directly from low to high momentum but must pass through intermediate 
scales giving a flow of particles from infrared to ultraviolet in our spectrum. 



Pi P2 




Figure 1: Elastic gluon interactions. 



Consider the elastic scattering of soft gluons, q± and q 2 , with hard gluons having 
Pi, p 2 on the order of T. The process is illustrated in Fig. 1 and gives a flow of 
energy from soft to hard as 



[2(iV 2 - l)] 2 f d 3 Pl d 3 p 2 d 3 qi d 3 q 2 



(2tt 



i 12 



^ 2E\ 2E 2 2uj\ 2uo 2 



{2nf5\p 1 + q 1 -p 2 -q 2 )-\M\ 



where 



with 



x [f q J pl (i + U)(i + f P2 ) - f q J P2 (i + f qi )(i + / P1 )] (wi - wa) 

'aN n 



M 



64tt 4 
M - 1 



7T 



2 r 



ut us ts 

o ~~ ~ ~~ 



t 2 



s = (pi + gi) 2 , t={q l - q 2 ) 2 , u = (p!- q x ) 



(7) 

(8) 
(9) 



The symbol R in the integration in (7) restricts u qi , u Q2 to be less than 
Now the hard particles p\ and p 2 will have a distribution given by (4) so 

1 + U = f P e E/T (10) 
allows one to write the bracket in (7) as 



[ ] = e E ^f Pl f p J q J q2 [(wi - u 2 )/T + 1/U ~ 1//, 



'/i J 



11^ 
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If f qi and f q2 also took the form (4), or (5), we would get exactly zero as expected in 
exact equilibrium. However, the incoming flux of soft gluons will naturally increase 
f q a little so that [ — uj 2 ) is greater than zero corresponding to a flow of energy 
from soft to hard particles. We can estimate that flow of energy by rewriting (7) as 



e 



_ [2(iV 2 -l)] 2 f _V?L ql ■ „ , (u) 



Taking — cj 2 to be positive (Eq. (12) is even under g x <-> g 2 ) and 



(Ecu) 



2 



Using 

from (2) and (3) and 



|M| 2 ~a 2 ^^- (13) 
1 1 m 4 v 1 

d 3 Pl ~ T 3 (14) 

d 3 gi~m 3 (15) 

g 2 ~ m 2 (16) 

c^i — cj 2 ~ m (17) 

along with (4) and (5) one gets 

e e <~^(— j^a/.lKim-^+r//, -%]-. (18) 

m/m ~ (amot) 1 / 4 (19) 
a/ g ~ aT/m ~ (20) 

from (2) and (5) one finds 

^ „ !!i a 9/4 (mot) 5/4 ^ _ ^ + _ T/u] 1 (21) 

Thus so long as m t 3> a~ 9//5 , f q will be very close to the equilibrium distribution 
making the [ ] in (21) small and compensating a 9 / 4 (m t) 5//4 . 

Although energy is transferred mostly in a direct way from the incoming particles 
of the source to particles having p ~ T there is also interaction with particles having 
momentum k where m <C k <C T. We can evaluate the rate of transfer of energy 
from (11) and (12) by the replacement P\,p<i — > &i,&2 along with the replacement 
d 3 pi — > d 3 ki ~ A; 3 in (14). One easily finds 

# ~ (™) f^-) 5 [ui-ot + T/U-T/UWT (22) 

V Ct. I lb J \ ' ' / 



after using (6). is the flow of energy from the particles coming from the source 
to those particles having momentum on the order of k. As we have seen earlier with 
the convention uj\ > uj 2 , [ } is positive so there is a positive flow of energy from 
the source to fc-particles, but this flow is suppressed by a factor k/T compared to 
the flow to particles of momentum T. However, since fk changes only very slowly 
this flow of energy must correspond to a flow of particles toward higher momentum. 
We can characterize this flow by a parameter, kfi ow , giving the average change in a 
particle's momentum with time in terms of which 

# = f k Ank% low . (23) 

In turn we can view ef as the flow of particles past a given momentum k or 

4 = kp k (24) 

where p k is the number of particles per unit volume which pass momentum k, from 
lower toward higher momenta, in a unit of time 

fjy 2 - 1) 

Pk = 2 (2^)3 f^k 2 k flow . (25) 

Using (24) and (22) one sees that p\ is independent of k so that there is a constant, 
in k, flow of particles from lower to higher momenta. One can check, using (22) and 
(24) along with (25) and [(wi - cu 2 ) + T/f q2 - T/f qi ] = m/(« 9 / 4 (m t) 5 / 4 ) from (21), 
that pk ~ dT 3 /dt so that the flow of particles is exactly the right amount to furnish 
the increase of the total number of particles of the system. 





P3 



Figure 2: 3 <-> 2 gluon processes. 

While the number of gluons in the system is growing with time p\ is much less 
than the rate at which gluons enter the system, — . Thus inelastic reactions are 
necessary in order to have chemical equilibrium. The dominant inelastic processes 
are those illustrated in Fig. 2 where gluons 1,2,3,4 have a momentum on the order 
of T and k <T. Dimensionally these graphs give (see Appendix A) 



T 6 ua 



3 



? nel ~ l-fpsfpti 1 + + /»)(! + fk) + UJvJkil + U(l + U\ (26) 
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or, using (2), (3) and (4) 



■inel „ 1^9/4^5/4 [(J _ T//J . I (2?) 

which is, parametrically, the same as (21). Again, when m < to < m one needs fk 
somewhat greater than T/u in order to transfer energy and lower the particle density 
to compensate for the gluons emerging from the source. We note that the LPM effect 
is small in the present case since the formation time of a gluon of momentum k is 

~ IF ~ (28) 

or 

with A the mean free path for the particle to undergo a collision of momentum 
transfer m, 

(30) 



1 TV 



A m 2 

Finally we evaluate the mean free path, A&, for a particle of momentum k to 
undergo a collision of momentum transfer k. Then, using (2), 

Ait 1 f k\ r 7 5-1 -i/4 



[a 7 (rnot) 5 ]- 1/4 , (31) 



t tk*fla 2 /k 2 \T 

so that we can expect the equilibrium spectrum fk — T/k for u > rh so long as 
m t > cT 1 ^ . The case m < u < rh will be discussed in the next section. 



2.2 The spectrum for (1/a) 7 / 5 < m t < (1/a) 9 / 5 

A major change occurs when m t = (l/a) 9//5 . We have seen that when m t > 
(1/a) 9 / 5 the gluon system has a equilibrium solution although there is a constant, 
in momentum, flow of particles from the infrared region to the ultraviolet region. 
However, from (21) it is clear that when m t < (l/a) 9 / 5 it is not possible to transfer 
the energy from the source in the region m < u> < rh to the momentum region near T 
fast enough using a near equilibrium distribution in the soft region. What happens 
is physically clear: The gluons from the source "pile up" in the region m < uj < rh 
until the occupation number becomes large enough to speed up the rate of transfer 
of energy from soft to hard modes so that the transfer can exactly compensate the 
rate of energy coming in from the source. For uj > rh in the distribution of gluons 



9 



3 




m 



in 



Figure 3: Gluon distribution as a function of frequency. 



will be very close to the equilibrium distribution (4), however, when m < uj < fn 
the distribution will be significantly altered. We now turn to a description of this 
altered distribution, illustrated in Fig. 3. 

We suppose that the source creates gluons uniformly in phase space in the region 
m < uj < in, with a constant inflow of energy e = m\ja defining the parameter m . 
This corresponds to an increase in occupation number f k ~ ^(^f ) 5 - Gluons having 
higher momentum, but located in the region m < uj < fn will elastically scatter 
with gluons having momentum on the order of T and will lose energy to the harder 
gluons as given by (18). Thus a gluon at point 1 in Fig. 3 will move to point 3 
further increasing f kz while f kl will be determined by the competition between the 
incoming gluons causing f kl to increase and the scattering with hard gluons causing 
f kl to decrease. Thus f k will increase rapidly as k decreases until f k becomes large 
enough that gluons having momentum k are absorbed by hard gluons, according to 
the graphs in Fig. 2, at the same rate that are created from the external source and 
from gluons which arrive at k after having elastically scattered with hard gluons. 
We now give estimates of these rates. 

Suppose &3 (point 3 in Fig. 3) is the momentum at which the loss by inelastic 
absorption, the graphs in Fig. 2, exactly balances the rate of arrival of gluons coming 
directly from the source or coming from the source via higher momentum regions. 
Then the rate of arrival of gluons is 



where — • — is the rate of the total number of gluons arriving over the whole 

a m, o o 



phase space m < uo < m and ( ^ ) expresses the fact that these gluons end up in 




(32) 
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the limited region of phase space on the order of fef. The rate at which gluons are 
absorbed by hard gluons, according to the graphs in Fig. 2 is 



f; 



abs 
k 3 



6^,3 



m J 



[fpJMi + U(i + U - /p 3 /p 4 (i + /pj(i + U(i + A,)] (33) 



or, using thermal distributions for f pi , f P2 , f P3 , f i 



PA 



fg 9 ~ -m^ f k , [- - l/f k ,\ T/m. 
Dropping the 1/ fk 3 term in (34) gives 

where we have anticipated k 3 /m <C 1 giving uj ~ m. 



(34) 



(35) 



Now we must also require that the rate of /c3-gluons cascading to even smaller 
momenta is not large compared to the rates given by (32) and (35). Suppose ak 3 - 
gluon goes to a momentum A; 4 via elastic scattering as given by Fig. 1, with the 
replacement qi — > k 3 ,q 2 — > k^. Then (see Appendix B) 



fi 



3^4 



-T 



a 



k,\ 



m 



ka m 2 \m J k 3 

which gives 

p3->4 



[UMl + U(l + f k4 ) - f P J ki (l + f Pl )(l + /*,)] (36) 



k 



^3 - ^4 

T 



+ l// fe4 - 1//; 



T 

-• (37) 
m 



Dropping the l/f ks and l/f k4 terms and using uj 3 - uj a ~ - ~ ^, since 

hi < k 3 <C to, one gets 



3^4 
k 3 



-my/af k3 ^/af k4 — — 



A:,i 



m 



m 



(38) 



Since we are looking for the value of k 3 at which the cascading to lower momenta 
ceases it is natural to require that (32), (35) and (38) be of the same size when fc 4 
is of the same size as k 3 . This gives, from (19), (32), (35) and (38) 



k 3 /m ~ [a 9/5 m t] 5/4 
V®fk 3 ~ [a 9/5 m t] " 5 



(39) 
(40) 



Finally, we give a calculation which shows that although f k increases sharply 
with decreasing k, in the region m < uj < to, it does not undergo a large change 
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as oj goes from above m to below m. To that end we first argue that in the region 
oj > m the occupancy f k lies not far from the thermal curve. Consider transitions 
^5 ~^ &3 by elastic scattering where 0J5 is just above m. Then, neglecting l/f ka 
compared to 1/ f ks (see Appendix B), 



5^3 
fe 5 



-m 



fe 5 



V®fk 3 — 



rn 



oj 5 — m 



T/m. 



The term v^/fe 3 (^) 3 is very large so that f k ^* 3 can be small only if 



(41) 



fk 5 ^ 



(42) 



oj 5 — m 

which is close to the thermal curve. Now consider transitions, due to elastic scat- 
tering on hard gluons, where k\ goes to k 3 . One has (see Appendix B) 



f, 



l->3 
k 3 



m 



Vaf kl Vaf k3 



m 



T 



T/m. 



(43) 



But the /fc 3 ~^ 3 in (43) cannot dominate (35) which means y/af kl cannot be large. 
Since (42) means \^af k5 is of order one, we see that there can be no strong change 
in fk as oj crosses the value m, despite the fact that the source turns on abruptly 
for oj < m. 



2.3 The spectrum for m t < (l/a) 7/5 . 

From (31) one sees that the mean free path divided by time for hard particles 
becomes less than one when mot < (1/a) 7 / 5 and thus that the majority of the 
gluons, and most of the energy, of the system are not in thermal equlibrium in the 
early time domain. Suppose, at a time t, the gluon spectrum has reached momentum 
Po(t). Then energy conservation gives 

m 5 



a 

or, using 



-t ~ f P0 pt> (44) 

m 2 ~ af po pl, (45) 
one gets 

p m ~ m\\/m§t. (46) 
For k/m^> 1 one expects f k to have the form 

r c(t)m 

fk ^ ■ 47 

aoj 
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The \jk dependence of fk will shortly be checked, but it naturally follows from the 
fact that the dominant interactions of gluons having momentum k are with gluons 
having momentum on the order of m and with momentum transfer of the order of 
m which is much less than k. This means the Boltzmann equation will obey 

y k (t)~m*Vlf k (t) (48) 

whose steady state solution is of the form (47). If we suppose that c{t) has a power 
dependence on mot, 

c(t) ~ a a (m tf, (49) 

then requiring fk(t) ~ l/a when m t ~ 1 with k ~ m and requiring fk(t) ~ 1/\fa 
when m t ~ oT 1 ^ and k ~ m, as given by (2) and (6), gives a = and b = —5/14 
so that 

c(f) ~ (m t)- 5/14 . (50) 

Eqs. (45), (47) and (49) give 

c(t) ~ (m t)- 5 / 14 (51) 



m 



while (46) and (50) give 



m ~m (m t) 1/14 , (52) 
p ~ m (m t) 3/7 . 



It is now easy to check that the mean free path, A Po , for gluons having mommentum 
p obeys 

^~ W^V^r 1 = ^ ~ 1 (53) 

as one might expect. We can now see somewhat more sharply why (50) must hold. 
If c(t) were to be parametrically larger than that given by (50) then X po /t would be 
less than one and the distribution (47) would not be stable. If c(t) were to be much 
smaller than given by (50) then the maximum value of fk which we would find, for 
very small values of k, would significantly exceed that given in (67) and would lead 
to the right hand side of (68) being much bigger than one, which is not consistent. 

Let us examine in a little more detail the consistency and stability of (47). In 
particular consider the emission or absorption of a gluon k, where m k po, by 
hard gluons as illustrated in Fig. 2. One has, as in (26), 

fk - [f P J P Jk(l + /ps)(1 + U - /ps/ P4 (1 + /*)(! + /»)(! + fk)] ■ (54) 
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Now the fs are all large so we can limit ourselves to terms with four factors of / 
on the right-hand side of (54). This yields 



fk ~ ~~ m?k? ^ k ^ Pl ^ P2 ^ PA f' P4 ^ ~ fpzfpiifpi ~^ fp^)] ~ f p\f pif pzf p^\ (55) 



Write 



E A = E 2 + u(l - A) (56) 



so that 



A ~ ^{(i - m 2 M^h + O + + fj} (57) 

where f = etc. In (54) we have done a dimensional analysis. The p^ factor 

is a rough way of writing the phase space factors coming from, say, d 3 pid 3 p 2 . Thus 
the fpiifp^fk + fp 2 ) in (57) should be understood as averaged over pi and p 2 . Since 
Pi and p 2 come in symmetrically we can rewrite (57) more simply as 

+ (58) 

where pi,P2 are on the order of po- Now the two terms, f'oufk and are of the 
same size if /& has al/k behavior as indicated in (47). Other power laws for fk will 
not allow fk in (58) to be zero. To see the stability of fk — suppose f' V2 ojfk + fp 2 
were positive. Then f k would grow in time. But this would increase f' P2 ujfk until it 
cancels f p2 . Thus (47) is a stable solution with respect to emission and absorption of 
gluons of momentum k by harder gluons. The exact value of c, in (47), is determined 
by the vanishing of < f' p uf k + fp > ■ 

Now let's investigate the stability with respect to elastic scattering k± + p\ — > 
k 2 + p 2 as illustrated in Fig. 1, and where m <C ki, k 2 <C p . One has 

fir 2 ~ [fM 1 + /p 2 )(! + /tei) - fM 1 + + (59) 

where we include dimensional factors for the pj and k\ phase space integrations, 
although we keep these momenta fixed in the Bose-factors to better understand the 
stability of the solution (47). Calling u kl — uik 2 = Acu, using (47) and writing 

/pi - f P , - -f P A" (60) 

fk, - f kl * (61) 
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one can write (59) as 

h 2 — 2 ^jT— [JmUJu + IpA • (62) 

Thus fl^* 2 will be zero under exactly the same conditions on c(t), in (47), as we 
found for stability under emission and absorption as given in (58). 

Finally, we turn to the region m < to < m. The situation here is very much like 
previously considered in Sec. 2.2. The rate at which gluons enter the system and 
then cascade down to a momentum fc 3 is given as in (32) 

where now m is given by (52). Eq.(35) is now changed to 



1 m 



(64) 



which is easily obtained from (58) by dropping f£ 2 compared to f' p2 ujf k . Finally (38) 
remains unchanged 

/r 4 ---v^/ fe3 v^/ fe4 (^) 3 (^), (65) 

and where the explanation of the meaning of this term is exactly as in the previous 
section. Requiring that the f k3 terms in (63), (64) and (65) be of the same size gives 



. 2 / \ 5 

m \ I m 



and 



h/m ~ ( - ] ( — ) ~ (m t)- 5/14 (66) 
\PoJ V m o/ 

ft, ~ i (^)* H 5 ~ i(mof) 15 /" (67) 

which now replace (39) and (40). We note that (66) and (67) agree with (39) and 
(40) when m t = a~ 7 ^ 5 . 

Finally, we note that the process where three soft gluons, having momenta on 
the order of k% as given in (66), annihilate into two harder gluons as illustrated in 
Fig. 4 with <7i , (72 5 93 ~ ^3 gives an f^ 3 of the same order as (63), (64) and (65). 
Indeed higher order processes where n gluons having momenta on the order of 
annihilate into 2 harder gluons also are of the same order (see Appendix C) since 

a (£Y/ fe ,~l. (68) 
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Q4 



Q5 



Figure 4: Inelastic 3 <-> 2 gluon interactions. 

3 Fixing the source in terms of wavelengths 

In the discussion of the previous section we have assumed that the source of incoming 
energy is spread uniformly in phase space in the frequency range m < uo < m where 
m and m are of the same order of magnitude. Now we are going to express the 
source in terms of wavelength rather than frequency. We again suppose that the 
rate of inflow of energy is 

e = ml/a, (69) 

but we here suppose that the gluons are uniformly deposited in phase space over 
the region 

< k < ho. (70) 

— * 

That is we take f source to be a constant, in both k and t, in the region (70) and 
where 

^t^f source J ° ^ k ^ k dk = m 5 /a . (71) 

In the present situation we have two parameters m and k with which to specify 
the source. (By allowing k to be time- dependent the current specification includes 
the source used in Sec. 2.) We shall, however, always suppose that k /m > 1 as 
the situation where ko is less than m seems to have an abnormally high rate of 
deposition of energy over a limited region of phase space. In the discussion that 
follows we deal with m and k as time-independent quantities, however, it will be 
easy to allow a slow time dependence of these quantities in our final results. 

Much of the discussion of the last section is unchanged since the basic pattern 
of the three different time regimes is set by the basic structure of the problem. We 
begin with the early time regime where mot < (1/a) 7 / 5 . From (52) and (66) we see 
that 

Y $ (™or 2/7 (72) 
k 
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and thus that k 3 /ko <C 1. Thus the picture here is essentially that described in 
Sec. 2.3, except that, when k /m <C 1, the region k Q < k < m now has no incoming 
gluons from the source. However gluons can come into this region from the process 
of four gluons, having momentum on the order of k 3 , going into a gluon, k, and 
two harder gluons. (We note that the three soft gluon to two harder gluon process 
cannot give a gluon, k, in the region ko < k <C m.) 

We can evaluate more precisely the dependence of when ko < k <C m, by 
requiring a cancellation in the rate of gluons of momentum k produced by the 
competing processes. First, the rate of production of gluons coming from qi + q 2 + 
<?3 + <?4 — ► k + k\ + kii with qi,q2,Q3, and q± of order k 3 while k\ and k 2 are of order 
m, is (see Appendix C) 

f£+ 3 ~mf k af m (73) 

where a (^-) 3 fk 3 ~ 1 has been used. Secondly the rate of gluons absorbed by hard 
gluons is, from (58), and after letting the 1/k 3 there become 1/m 3 reflecting the 
cutoff of the 1/k 3 singularity when k/m <C 1, 

ff- ~ mh t ( 1/A - ™) (74) 

where we have used 

< f p c(t)m/a + f p >= (75) 

the stablity condition determining c that we found in Sec. 2. 3. In arriving at (74) 
we have also used (fj^) 2 cxf Po ~ 1 and c = m/po- Finally, there is the rate at which 
gluons k convert to gluons of momentum around k 3 by elastic scattering with hard 
particles, 

ft" 3 ~ -P 3 oS ( y IfMl + fk 3 )(l + U ~ fk 3 f P2 (l + /*)(! + f Pl )} ■ (76) 

7TI \ 7TI J K 

It is straightforward to get 



which, after using (75), becomes 



a{U - m) -l/fk 



m c 



(78) 



If fk were near its normal size as given by (47), then (78) would dominate. (Re- 
call c/ ^fa ~ [a 7 / 5 m t] 5 ^ M is a large quantity in the time domain we are currrently 
considering.) This means that 1/ fk must be much smaller than "normal", and the 
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1/ fk term can be dropped in (74). The remaining term in (74) is comparable to 
(73) so that one can get stability by requiring 



—ua/m 



m 
~k 



a(cu — m) 



m c 



-1/fk 



or 



fk 



mc 



a u -m + O([m t]- 5 / 14 £:) 
so that fk begins to rise rapidly as k decreases beyond m. 



(79) 



Now turn to the region m t > (l/a) 7 ^ where gluons having w/m > 1 are in 
thermal equilibrium. From (39) 



k 3 /k ~ [a 9 ^m t\ 5 ^ m/k . 



(80) 



So long as k 3 /k <C 1, that is so long as m t <C (l/a) 9 ^(k /m) 4 ^, the large peak in 
fk will still occur when k ~ k% and (39) and (40) should remain valid. Let us again 
investigate fk in the region fc < k < m when k^/ko <C 1. Now the three relevant 
processes are 



fk~ 3 ~ rnfkafr, 



ah 



in 



mfkOifm [a 7/5 m t] 



while 



ft" 3 ~ ~h [a 7/5 m t] 



-5/4 



uo — m 



-l/fk 



T/m 



(81) 



(82) 



and the process where k is emitted or absorbed by thermal particle as illustrated in 
Fig. 2, 

(83) 



UJ 



ft ~ -mfkV^ - l/fk\ T/m . 

The dominant contributions now are (82) and (83). Requiring a cancellation to get 

a steady state gives 



fk ~ — [(w - m)/m + O(A;/m[a 9/5 m t] 5/4 )] 1 . 

7TZ 



(84) 



Thus /fe grows as k decreases below m. When mot approaches (l/a) 9 ^ 5 (ko/m)^ 5 
the picture starts to change as k 3 approaches k . At m t = (l/a) 9 / 5 (k /m) 4 / 5 the 
occupancy f k begins to have little k— dependence for k/k < 1. Requiring f* ource ) 
as given by (71), to cancel with f^ bs , as given by (83), in the domain < k < k 
gives 

(85) 



V^fk~ [a 9 / 5 m t]' 5/ \m/k ) 3 
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We note that (84), evaluated at k = k is much smaller than (85). Thus for m t > 
(l/a) 9 / 5 (k /m) 4 / 5 there is a "discontinuity" in f k with 

Vufk\k=h -e ~ V&fk\k=k +e ~ [a 9/5 m t] 5/4 (m//c ) 3 . (86) 

By the time mot ^> a~ 9 / 5 (m/ko) 12 ^ 5 This jump in fk at ko has become very small, 
compared to the size of fk, and complete "equilibration" is reached. Thus if ko is of 
the order of m then complete thermalization occurs at a time m t ~ (1/a) 9 / 5 as in 
Sec. 2. If ko/m <C 1 then complete thermalization occurs when 

mot ~ (l/a) 9/5 (m/A;o) 12/5 . (87) 
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A The derivation of Eq. (26) 

In this Appendix we consider the 2^3 gluon interaction where 4 gluons have 
momenta on the order of T and the momentum of the remaining gluon is k <C T. 
To derive the result given in Eq. (26), we follow Ref. [22] where the 2^3 process was 
calculated in all details. The most relevant process is the large angle gluon-gluon 
scattering with a small angle gluon radiation/absorbtion (gluon-gluon scattering 
with collinear emission/absorbtion). The inelastic 2 <-> 3 process is given by 

[2(iV c 2 - l)] 2 f d 3 Pl d 3 p 2 d 3 P3 d 3 P4 d 3 k , , 4 . 4 , 



x|M 1M4l | 2 [f r J r J t (\ + /„,)(! + /„) - /„/„(! + / P1 )(l + /„)(! + f k )} 



where N = 2(iV c 2 — 1) J -^sfk and the symbol R in the integration in (88) indicates 
that k Pi,P2,P3,P4- To evaluate (88) one supposes that the 5-gluon process 
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consists of a hard collision involving four gluons. Then the fifth gluon, k, is collinear 
with one of the incoming or outgoing gluons, giving a factor of 4 which when devided 
by the symmetry factor, 2, explains the appearance of 2 in (88). Since all four 
possibilities lead to the same result, let us consider only the case where the gluons 
4 and k are collinear. Collinear means that the invariant mass = (p 4 + k) 2 is 
significantly less than Q 2 = min(|s|, \t\, \u\), where Q is the momentum scale of the 
hard scattering. 

In Ref. [22] (see Eq. (20)) it has been shown that the following factorisation holds 
in the collinear limit 

f d 3 p 3 d 3 p 4 d 3 k 4 4 2 

J WW WW WW (27r) 5 fa + - k ) \ M ^\ 

f d 3 p 3 d 3 p 4 4 4 2 f dzds 4k a s (s 4k ) 

^JWWWW {2n) 6 i m «-«i x Jw k — *r PM 

(89) 

where P gg (z) is the unregularized g — > gg splitting function 

P 99 (z) = 2N C + _iL_ + *(i - z )) ( 9 °) 

and |M 12 ^34| 2 is the matrix element for the hard scattering given by (8). With (89) 
Eq. (88) becomes 

[2(iV c 2 - l)] 2 f d 3 Pl d 3 p 2 d 3 p 3 d 3 p 4 , o . 4 _ 4 , 



N = 2 W^ J R 2f 2f 2f 2f (27r) 4 <Hpi+P2-P3-P4) 

X / — [/p 1 /p 2 / fe (l + /p3)(l + /pJ-/p 3 / P4 (l + /p 1 )(l + /p 2 )(l + / fe )]- 

S 4 fc 27T 

(91) 

The dominant contribution to |Mi2^3 4 | comes from the small-t region 

64tt 4 (aN^ 2 



W - tt-V— - I — ) -y 02) 



where « ~ -s for small £ = {p\— p 3 ) 2 due to s + w + t = 0. Using also s = (pi+p 2 y 
2E 2 , the integration over p 4 in (91) gives 

N = m l~^ )]2 (—) 2 [ d 3 Pl d 3 p 2 d 3 p 3 5\E, + E 2 -E 3 - E 4 ) 



4(2tt)« V * J Jr ------- v i( Pl - P3 ) 2 } 2 

x [ < ^n ( ^Ap gg{z ) [/ pl / P2 / fc (i+/ P3 )(i+/p 4 )-/p3/p 4 (i+/p 1 )(i+/p 2 )(i+/ fe )]- 

J R S 4 fc Z7T 

(93) 
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With d 3 p 3 = d 2 p 3 j_dp3L = d 2 p 3 j_dE 3 , where p 3 ± is the transverse and p 3 i the longin- 
tudinal component of p 3 , and (pi — p 3 ) 2 = p 2 ± , one can perform also the integration 
over dE 3 and obtains 

■ [2(N 2 -1)] 2 {aN c V f j3 j3 j2 1 

dzds^ajp^ [U n f k (l+S n ){\+f„)-f r J n {-i+f n )0-+U^+h)\- 
R s 4fc Ztt 

(94) 

Now, with pl ± ~ m 2 because of screening, d 2 p 3 ± ~ m 2 , d 3 pi ~ T 3 , d 3 p 2 ~ T 3 and 
the relation e = Nw, one obtains the parametrical estimate given in Eq. (26). Note 
that k <C Pi,P2,P3,P4 is satisfied by requiring z to take only values close to in (94). 



x 



B The derivation of Eqs. (36), (41) and (43) 

In this Appendix we show how to get the parametrical estimates given in Eqs. (36), 
(41) and (43). Let us first derive Eq. (36). The elastic scattering of soft gluons, k 3 
and ki {k 3 > k^), with hard gluons having momenta p±, p 2 on the order of T, as 
shown in Fig. 1 with the replacement qi — > k 3 and q 2 — > k 4 , is described by 

. 4 2(iV c 2 - 1) f d 3 Pl d 3 p 2 d 3 k 4 4 2 

x [/J B (1 + /*.)(! + U - fU P2 (l + / fe3 )(l + • (95) 

where \M\ 2 is given in Eq. (8). The hard momenta read pi = (Ei,pi) with Ei m \pi\ 
while the soft momenta are kj = (cuj, kj) with oj 2 = m 2 + k 2 ~ m 2 . 

It is the small angle scattering (t small) which gives the dominant contribution 

to |M| 2 , 

,,, l2 64tt 4 (aN<> 2 
\M 2 - 



MSI , , 

(A? - 1) V W it!' (96) 

Taking into account that t ~ m because of due to s + -u + t = 0, 
and s — (pi + k 3 ) 2 « 2Eio;3, one obtains 

IMI 2 ~ ^ (" N 'X {2Em? (97) 

Using also 

S\pi + h-p 2 - fc 4 ) = 5 3 (p! + k 3 -p 2 - h) 5(E 1 +uj 3 -E 2 - uj 4 ) (98) 
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the integration over p 2 in Eq. (95) then gives 

Jk 3 ~ 



... .., , tffc4 +w 3 ~y(Pi + k 3 -fa) 2 -u A ) 

x[/ fc3 / P1 (l + / fe J(l + / P2 )-/ fc4 / P2 (l + / fc 3)(l + / P1 )] , (99) 



where / P2 within the brackets [ ] is to be evaluated at p 2 = pi + fa — fc 4 . 
The evaluation of the delta-function 



6(E 1 + u 3 - \J (pi + fa - fa) 2 - uj 4 ) 



8(\pi\ - m + — - - |pi| - -r=n(fa -fa)-m- — 
2m pi 2m 



5(-\fa\cos(9)) , 



(100) 



where 6 is the angle between pi and fc 3 , allows us to perform also the angular 
integration, d 3 pi5(— \k 3 \ cos(0)) = 2np 2 d\pi\/\k 3 \, which leads to 



3^4 

k 3 



16 (aN r ) 2 f 2 , 3 , 1 
Pid\pi\ a fa —r- 



(2tt) 2 m± J - — ' \ ksl 
x [fkJ Pl (l + /*.)(! + U - fkJ P2 (l + / fc3 )(l + f Pl )\ ■ (101) 



Now, with p\ ~ T 2 and d\p\\ ~ T, is easy to see that Eq. (101) is parametrically 
the same as Eq. (36). 

To get Eq. (41), one follows the same procedure as above, with the replacement 
fa — > fa and A; 4 — > fc 3 , and gets 



^5^3 
■U 5 



16 (aiV c ) 2 f 2j _. j3; 1 



(2tt) 2 m 4 J — " ' - |fc 5| 
x [fkJ Pl (l + /*,)(! + /«) - + /*»)(! + /*)] 



Using (10) which allows us to write the bracket in (102) as 

[ ] = e E ^ T f Pl f P2 f k Jk 3 [(w 5 - cu 3 )/T + l// fe3 - l// feB ] 
and pi,p2 ~ T, <i 3 A;3 ~ fcf and (3), we get 

3" 



(102) 



(103) 



5-*3 



~ — m 



V®fk 3 — 



in 



u 5 - u 3 



T/m (104) 



which for m/fa ~ 1, u> 3 ~ m and l//fc 3 <C l//fc 6 reduces to the result in Eq. (41). 
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To get Eq. (43), one starts with 



l-t-3 
k 3 



2(n 2 c - 1) r d 3 Pl d 3 P2 d 3 h 

(2tt) 5 2lu 3 J ~2E[ ~2E 2 2^7 



S 4 (pi + k 1 -p 2 -h) \M\ 



x + / fe3 )(i + U - f k J P2 (i + AJ(1 + /*)] 

and follows the same steps as in above calculations which lead to 

lj i — m 



l->3 
fc 3 



my/af kl V<xfi 



k 3 



T 



~l/f> 



fci 



T/m , 



(105) 



(106) 



where the factor ik\jm) 2 comes from the phase space d 3 k\ and the angular inte- 
gration d 3 p\ S(— \ki\ cos(6>)) = 27rp| d|p2|/|^i|- Since k\ is close to m in the case 
discussed, it is legitimate to set (ki/m) 2 ~ 1 which gives the estimate in Eq. (43). 



C n —> 2 gluon merging 

Here we show that the process where n soft gluons, having momenta on the order 
of ks as given in Eq. (66), merge into two harder gluons ( with ki, k 2 ~ m) gives 
an fk 3 on the same order as (63), (64) and (65). It is instructive to begin with a 
simpler process where only 3 three soft gluons annihilate into two harder gluons as 
shown in Fig. 4, with qi,q 2 ,q3 ~ k 3 , and then go over to the more general case. For 
the 3^2 process we have 

3^2 2(N 2 - 1) f d 3 Ql d 3 q 2 d 3 h d 3 k 2 ,4/ , , ; / \ I ,,|2 



' fc3 (27r) 8 2^3 7 2 Wl 2a7 2 2cu fcl 2 W)fc2 W 2; 1 

x [/JJJl + + / fe2 ) - / fel / fe2 (l + + U)(l + /»)] (107) 

and after the dimensional estimate 

|M 2 |~a 3 A7 (108) 

and the integration over A; 2 

fl^ 2 ^ ~f J rf3 ^i rf3 ^2 rf3 ^i S(u qi + Wg 2 + - id kl - t^ 2 ) 

x lUfMl + f kl )(l + fk 2 ) - f kl f k2 (l + + U(l + f q3 )} (109) 

Doing another dimensional estimate 

d 3 kx 5{id qi + id q2 + id qz - id kl - id k2 ) ~ m 2 , (110) 
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keeping only the leading contribution in the brackets 

lUfMi + + fk 2 ) - fMi + U(i + U)(i + U] 

= [fks fm] + corrections 
and using <i 3 gi, d 3 q2 ~ fcf, we find 



fir 2 ~ m f m afk 



a 



fk 3 



which with (68) reduces to 



(111) 

(112) 
(113) 



and is of the same order as (63), (64) and (65). 

Following the same procedure, it is easy to show for the n — > 2 process that 



n->2 
k 3 



mf m af k 



a 



fk 3 



(114) 



which again, when using (68), becomes of the same order as (63), (64) and (65). 
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